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Abstract 

A model function relating effective stress to fracture permeability is developed from Hooke's law, 
implemented in the tensorial form of Darcy's law, and used to evaluate discharge rates and pressure 
distributions at regional scales. The model takes into account elastic and statistical fracture parameters, 
and is able to simulate real stress-dependent permeabilities from laboratory to field studies. This modeling 
approach gains in phenomenology in comparison to the classical ones because the permeability tensors may 
vary in both strength and principal directions according to effective stresses. Moreover this method allows 
evaluation of the fracture porosity changes, which are then translated into consolidation of the medium. 
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Introduction 

Crystalline and non karstic sedimentary rocks are 
anisotropic geological media with low hydraulic con- 
ductivity |22j . In such media, groundwater flow oc- 
curs primarily and sometimes exclusively through 
non-filled fractures. Their spatial arrangement (i.e. 
fracture network) leads to groundwater flow at a re- 
gional scale. At this scale, the most simple and useful 
way to conceptualize these aquifers is the equivalent 
porous media; the principal permeabilities of each 
fracture family are combined in space and result in 
a tensor describing the equivalent hydraulic conduc- 
tivity of the rock mass |13| [2] . 

The sensitivity of aquifers dynamics to effective 
stress was first described for granular porous media 
[26]; the process was then also observed in fractured 
aquifers [l9]. Nowadays, the dependency of fracture 
permeability on effective stress is a well known re- 
search topic and has been intensively studied during 
the last decades, especially to evaluate the stability 
of rock masses in presence of dams, tunnels, geologic 
radioactive waste repositories or C02 sequestration 
fields pSl [m [231 inni E] . in regional and deep ground- 
water flow systems, the reduction of water pressures 
leads to increasing effective stresses and decreasing 
permeabilities, with a possible consolidation of the 
aquifer. On the contrary, increasing groundwater 



pressures result in decreasing effective stresses and 
in increased permeabilities. Based on field and lab- 
oratory test results, Louis |19j and later Walsh [28] 
derived respectively an exponential and a logarithmic 
model to explain permeability decreases with increas- 
ing effective stresses. The relationship between effec- 
tive stress and permeability has been clearly identi- 
fied both at local and regional scales via laboratory 
tests (see [271 E E] ) , field tests (see [1 [25]), and ob- 
servations of aquifer consolidation by measurements 
and modeling of ground subsidence (see [HI [231 [30] ) • 

However, the equivalent porous medium ap- 
proach and consequently the classical Darcy solu- 
tion implemented in regional groundwater numeri- 
cal models generally ignores this relationship [21j . 
hence the interest in modifying the flow equation 
to a more realistic one, explicitly accounting for 
stress-dependent permeabilities. The present ap- 
proach consists in inserting constitutive laws relat- 
ing effective stress to permeability in the tensor form 
of Darcy's law, so that the permeabilities vary with 
stress (depth and geology) and water pressure. The 
constitutive model must: (1) respect most of the 
physical process at the microscopic scale but should 
also lend itself to practical application at large scales, 
and (2) be simple from a numerical point of view. 
Note that a number of rock mechanics codes exist 
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(see for example [TH [H [5l [29] ) that solve coupled hy- 
dromechanical problems. However, these generally 
apply to relatively small scale problems, because they 
involve full and detailed deformation processes, and, 
therefore, become computationally prohibitive at hy- 
drogeological scales. On the contrary, the present 
work focuses on a macroscopic approach allowing ef- 
ficient large scale computations, while preserving the 
essence of the hydromechanical processes. 

Firstly, a constitutive model is presented where 
fracture permeability is a function of the effective 
stress, as well as of the statistical distribution of 
the length of the asperities and their elasticity. Ex- 
pressed in its tensor form, this law describes the pro- 
cess at the rock mass scale. Secondly, simulated 
stress-dependent permeabilities are compared with 
laboratory and field measurements from Durham [7] 
and Cappa [Ij. Thirdly, a finite element simulation 
is performed in order to illustrate this modeling ap- 
proach. The constitutive model is also used to eval- 
uate the changes in porosity between an initial and a 
modified hydrogeological state, and to compute the 
resulting subsidence. 

Constitutive aperture-stress model 

The model considers a single fracture as a pair of sur- 
faces, characterized by a set of asperities, the length 
of which follows a statistical distribution. This asper- 
ity population can be characterized by fracture mor- 
phology analysis (see O HO]). Assuming that each 
asperity i obeys Hooke's law, the resulting normal 
stress, (Tj, proportional to its deformation is: 



Azi 



Ei 



(1) 



where the symbols stand for asperity original length 
Zi, compression Azj = — a, elastic modulus Ei, av- 
erage asperity section Sj, exerting force Fi and frac- 
ture aperture a. Eq. ([T]) implies the following condi- 
tions: 

First, if a > Zj, (Tj = (the asperity is at its original 
length) . 

Second, if a = 0, fjj = Ei (the asperity is subjected 
to a total compression). 

Under a given normal stress, the asperities return 
an equilibrium equivalent stress a defining a specific 
aperture a (Figure [T^). This normal stress results 
from the integration of all stresses exerted by indi- 
vidual asperities compressed to various degrees. For 



a given aperture, the probability that an asperity is 
in contact with both fracture faces P{zi > a) cor- 
responds to the ratio between the number of com- 
pressed asperities and the total number of asperities: 



P{zi > a) 



D{z)dz 



Nt 



(2) 



where D{z) is the statistical distribution of the as- 
perity lengths (Figure [Tja), is the number of com- 
pressed asperities and Nt is the total number of as- 
perities. In Eq. ([2]) the infinite upper bound of the 
integral can be replaced by the maximum fracture 
aperture oq, which also represents the original length 
of the longest asperities. Note that the integral of 
D{z) must be equal to unity. Glover et al. [10] note 
that D{z) is frequently assumed of Gaussian type. In 
this paper, a number of simple typical distributions 
are considered, as well as the more realistic Weibull 
distribution. 

Assuming average values for asperity elastic mod- 
ulus and section, and associating [2] and [l| the equi- 
librium normal stress for a fracture with an aperture 
a is obtained by weighting each asperity contribution 
by its probability density. Integrating over all active 
asperities yields: 
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where F is the force exerted by the compressed as- 
perities, A is the fracture surface area, E is the elas- 
tic modulus of the fractured rock and s is an aver- 
age asperity section. Nt/A = r] is the asperity areal 
density. Eq. Q respects the same conditions as 
Eq. Firstly, the maximum fracture aperture 

ao is reached at cr = (a = ao, no compression). 
Secondly, total fracture compression (a = 0) occurs 
when a = ijEs = do, where cto is the fracture closure 
normal stress. 

Model adjustment for different statistical 
distributions 

Different cr(a) models are obtained depending on the 
statistical distribution D{z). For example, for the 
uniform distribution (Figure [2^) D{z) = l/ao , Eq. 
^ becomes: 
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Figure 1: (a) Schematic illustration of a rock mass intersected by a fracture and its conceptualization with a set of 
asperities. The fracture under the normal stress a has the aperture a, the maximum fracture aperture oq is reached 
when there is no stress, (b) Continuous statistical distribution, D{z), of the asperity length, z, and probability of 
contact. 



Several constitutive models are found proceeding in 
the same way for different types of distribution D(z) 
(Table 1). After integration the fracture aperture a 
is directly related to normal stress a. As mentioned 
in the introduction, the model must be simple from 
a numerical point of view, hence the need to refor- 
mulate the specific models presented in Table 1 in a 
generic equation of the form: 



ao 



,n> 1 



(5) 



Depending on the value of coefficient n, Eq. ^ pro- 
vides exact stress-dependent apertures for the non- 
logarithmic functions of Table 1 and good approx- 
imations for the logarithmic ones (Figure 2b). The 
symbol n stands for the coefficient of asperities length 
statistical distribution. Statistical distributions char- 
acterized by many large asperities, such as singular 
and linear increasing, get low coefficients n, 1 and 
2, respectively. On the contrary, distributions with 
many small asperities (linear decreasing and Weibull) 
are correctly approximated with relatively high co- 
efficients n (4.7 and 9). Inversely, the distribution 
D{z) can be found for a given stress/aperture func- 
tion a{a). Differentiating Eq. (3) twice with respect 
to a yields: 

d'^a{a) (To 
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For the general model assumed in Eq. ([s]), this yields: 
d'^a{a) n{n — 1) 



da'^ an 



n-2 



(7) 



and the distribution D{z) is obtained by equating 
Eqs. ([6]) and ([7]): 



D{z) =n(n-l)4^ ( 1- — 



n-2 



,n>l (8) 



Relation with hydrogeological parameters 

Because of the saturated flow conditions considered 
in this work (i.e. fractures are completely filled by 
water exerting the pressure p), the normal effective 
stress a' is taken into account instead of the normal 
total stress a. In the case of lithostatic stress condi- 
tions az = Pt-gZ and in the absence of shear stresses, 
the resulting effective stress a' acting perpendicularly 
at a depth Z on a given fracture plane is obtained by: 



a = crn • n — p 

,2 , \^2 



a^X 
cr^A 
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where pr is the rock mass density, g is the gravita- 
tional acceleration, nx,ny,nz are the components of 
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Table 1: Example of aperture-stress models for different statistical distributions of asperity length. 
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Figure 2: (a) Possible statistical continuous distributions of asperity length in a fracture, and (b) corresponding 
stress/aperture solutions, with their approximation by Eq. ([5]). 

the unit vector n normal to the fracture plane, pw ratio of horizontal to vertical stress, 
is the water density, h is the pressure head, and a 
is the Biot- Willis coefficient. The A coefficient is the 
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Expressing Eq. ([s]) for the aperture a: 



a = ao 



(10) 



and assuming the vahdity of the cubic law in the 
fractured rock, the stress-dependent permeabihty is: 
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:ii) 



yielding the hydrauhc conductivity parallel to frac- 
ture plane: 



where: 



a 



PwQ fal 



(12) 



(13) 



p-w 12 

with a maximum Kq for a' = 0. The symbol 
/ = Nf/d is the frequency of the fracture family, 
namely the number of fractures Nj counted over a 
distance d, and is water viscosity. Note that 
Eq. (12) is very similar to the constitutive models 



proposed by [171 [l5] . The same model function was 



found by [9] via a different approach. Eq. (12) can 
be used to compute the equivalent macroscopic hy- 
draulic conductivity tensor of a rock mass intersected 
by m fracture families using the tensor summation: 
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for each fracture family i, Kq. is the maximum paral- 
lel hydraulic conductivity, cr^ is the normal effective 
stress, (Tq. is the fracture closure normal stress, rij 
relates to the asperity distribution, / is the identity 
matrix, rij is the unit vector normal to the fracture 
family i, and (g) denotes a tensor product. 

If the contribution of the rock matrix is neglected, 
the porosity (p of the fractured rock mass is: 



i=l 



(15) 



Introducing Eq. (10) into Eq. (15), a stress 



dependent porosity is obtained: 



E^o. 

i=l 



cr, 



0, 



(16) 



Always neglecting the contribution of the rock ma- 



trix, Eq. (16) can be introduced in the definition of 



the specific storage coefficient, 5s: 
Pwg4> 



5, 
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where 0o and 
ity and maximum specific storage coefficient, respeo 



SQ are the maximum fracture poros- 



tively. Eq. (16) can be used to determine, for each 
fracture family i, the vertical variation in fracture 
porosity, Ai;^, due to a change in effective stress, 
between an initial and a successive hydrogeological 
state at elevation z: 



1=1 

m 

= E^o 



1=1 



0, 



(18) 



where the symbols (ph and (phs stand for fracture 
porosity at an initial and at a successive pressure 
head state. The multiplication with the component 
Uz of the unit normal vector n is used to obtain the 
vertical change in fracture porosity. Integrating all 
the porosity changes in the vertical direction, from 
the bottom of the aquifer Zj, to the top zt , results in 
the local settlement: 



T{x,y) 



A(f){z) dz 



(19) 



Eq. (19) provides: (1) where T > is aquifer verti- 
cal consolidation under increasing effective stress and 
(2) where T < is aquifer vertical expansion under 
decreasing effective stress. Finally, considering both 
the hydraulic conductivity and the specific storage 
coefficient as functions of effective stress results in 
the non-linear groundwater flow equation: 

Ss{a')— = V-{K{a')VH) ; H = h + z (20) 

where H is the hydraulic head, K((j') is the hydraulic 
conductivity tensor as expressed in Eq. (14), 5s((t') 



is the specific storage coefficient as expressed in Eq. 
(17), and t is time. The symbols h and z stand for 



the relative pressure and elevation head, respectively. 
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Figure 3: Verification of Eq. (12 1 model by comparison with stress-dependent fractured rock permeabilities of (a) 
Durham [T and (b) Cappa [1]. 



Comparison between simulated, experi- 
mental and field measured permeabilities 



Eq. (12) is verified by comparison with stress- 
dependent permeabilities from Durham [7] and 
Cappa [4j. Laboratory tests carried out by Durham 
[7] showed the behavior of the permeabihty of a frac- 
ture sample, taken at approximately 3.8 km depth, 
when subjected to an increasing confining pressure 
(stress). Simulated permeabilities correspond well 
to those measured by Durham [7j, especially for 
high stresses (Figure [3^, Appendix II). With exper- 
iments at shallow conditions, Cappa |4] investigated 
the pressure-dependent increase and decrease of frac- 
ture aperture. Results showed that fracture aper- 
ture is subjected to hysteresis process. Eq. (12) is 



used to fit the field data of Cappa [3] (Figure |3|; 
which for that example were transformed from aper- 
ture and water pressure to permeability and normal 
effective stress. Fitted parameters are given in Ap- 
pendix II. Also for this example, the model provides 
a good comparison between simulated and measured 
data. However, only the rising branch of the hys- 
teresis curve, which corresponds to an increasing wa- 
ter pressure and a decreasing effective normal stress, 
is correctly simulated. Eq. (12) can not reproduce 
a hysteresis, because the model describes only the 
elastic part of deformation, and therefore simulated 
permeabilities will be the same for rising or falling 
effective normal stresses. In Figure [3|3 the hysteresis 
occurs because the tested rock does not exactly follow 



Hooke's law. However the variation of permeability 
is so low, that this phenomenon may be neglected at 
regional scale. Note that, numerical values of the co- 
efficient of asperities length statistical distribution, 
n, and the fracture closure effective stress, (Tq, can 
be obtained by calibration of Eq. (12) on measured 



stress-dependent permeability data, and applied for 
large scale analysis. 

Illustrative examples 
Steady state 

The preceding equations were implemented in the 
multi-purpose Groundwater finite element software 
[6], in order to illustrate: (1) how simulation re- 
sults vary, if the effect of effective stress on hydro- 
geological parameters is taken into account; (2) the 
regional effects of a deep tunnel. At steady state 
and with stress-dependent hydraulic conductivity the 
flow equation is: 

V • (K((j')Vif) = ; H = h + z (21) 

The virtual model domain is a 2D vertical cross sec- 
tion representing an Alpine hydrogeological system 
composed of three geological formations, completely 
saturated with water, with different hydraulic prop- 
erties (Figure [4| Appendix II). Note that, in Ap- 
pendix II are shown numerical values used in sim- 
ulations, these are based on field investigations of 
the Emosson fractured rock mass (Switzerland), us- 
ing the method described by Kiraly [11] . In appendix 
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II, Kmax and Kmin are the eigenvalues of the hy- 
drauhc conductivity tensor, and 9 is the angle be- 
tween the horizontal plane and the direction of K^ax- 
The stress field is defined by vertical stresses a{z) set 
equal to the lithostatic pressure (the weight of over- 
lying rocks above elevation z): 



a{z) = g Pr{u) du 



(22) 



In the present study horizontal stresses (Jx are 1.5 
times stronger than the vertical stresses. This ap- 
plies well to orogenic belts or areas that have been 
glaciated, such as the Alps (see |20]). Before tun- 
nel construction, a steady state flow is assumed from 
the highest points (crests) to the lowest points (val- 
leys), by specifying boundary conditions at the do- 
main surface as atmospheric pressure {H = z), and 
at other limits as no-flow conditions. Then, a tun- 
nel is constructed. A constant atmospheric pressure 
is specified in the tunnel indicating that it behaves 
as a draining structure and consequently increases 
the effective stress which causes aquifer consolida- 
tion. Several simulations are computed to compare 
the present approach with the classical one neglecting 
the dependence of permeability on effective stress, 
and to study the influence of the coefficient n on 
discharge rates, pressure and consolidation distribu- 
tions. Consolidation is computed between the initial 
state (without tunnel) and the disturbed state with 
perturbation caused by the tunnel (Figure [s]) . This 
illustrative model is directly inspired by real cases 
of fractured aquifer consolidation caused by tunnels 
excavation (see [16j, and [30j). 

Results and discussions 

Results show that the introduction of stress- 
dependent permeabilities in Darcy's law leads to 
lower discharge rates, relative to the classical ap- 
proach that only considers constant permeability, es- 
pecially for high values of coefficient n (Figure 6a). 
This reduction in discharge rates is directly related to 
the variations of the hydraulic conductivity tensor in 
both strength and principal directions according to 
effective stresses. This spatial variation of hydraulic 
conductivity tensors also impacts the distribution 
and the magnitude of hydraulic heads, flow paths, 
flow velocities and transit times (Figures [5]and|6]3). 
Overall permeabilities decrease in the deeper areas 
of the domain, while they tend toward the maximum 
near the surface. As previously mentioned, a high 



value of the coefficient n indicates a predominance of 
relatively small asperities. In such a case, the drop 
in permeability will be significant, because there are 
only a few large asperities to oppose the increase in 
normal effective stress. Compared to the classical ap- 
proach, the impact of the tunnel on the system ap- 
pears weaker when considering stress-dependent per- 
meabilities. For the proposed method, the highest 
consolidation occurs in systems with an intermedi- 
ate n value (1 < n < 5), because they are the most 
sensitive to pressure change with the largest poros- 
ity variation (Figure [6^) . Fractured systems featur- 
ing large asperities are less affected by the process, 
because the asperities stop the closure. Overall, the 
magnitude of aquifer consolidation is low because the 
proposed method computes only elastic reversible de- 
formations obeying Hooke's law, and acting on frac- 
ture network porosity. Moreover, boundary condi- 
tions specified at the domain surface provide unlim- 
ited water inflows that dampen aquifer depressuriza- 
tion. 

Transient state 



Eq. (20) is used to solve the transient groundwa- 



ter flow problem having the same model domain, 
boundary conditions at the tunnel and hydrological 
parameters. The initial hydraulic heads are taken 
from the steady state model without tunnel. On the 
domain surface two different types of boundary con- 
ditions are tested: (1) constant atmospheric pressure 
{H = z); (2) no- flow condition. This no- flow con- 
dition could represent an aquifer filled with connate 
pore waters and isolated from recharge zones, or a 
confined aquifer suddenly cut-off from its recharge 
zone. 



Results and discussions 

For the first case where temporally constant atmo- 
spheric pressure hydraulic heads are specified at the 
domain surface, the initial and final discharge rates 
as well as the vertical settlements match those simu- 
lated by the steady state models. In transient state, 
the tunnel causes a hydraulic depressurization of the 
rock mass followed by gradual aquifer consolidation 
(Figure [7^). For the second case with sudden no-flow 
condition at the domain surface, the tunnel drainage 
empties the system, which becomes hydrostatic. The 
recession curve of the water drained by the tunnel 
rapidly runs dry. In such a case, the magnitude of 
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Figure 4: Model domain and boundary conditions; the hydrogeological system is composed of three rocks, each one 
exhibiting a different fracture network. Rock 2 is the most permeable, while rock 3 is the least. Three simulations are 
run for each statistical distribution of the asperities length: one without tunnel (natural system) , one with tunnel, and 
finally one that simulates the aquifer consolidation. 




Figure 5: Hydraulic head, flow paths and infiltration/exfiltration fields for the classical approach, and for the Weibull 
distribution at initial conditions (a, b) and after the tunnel introduction (c, d). Note that, in (b) and (d) fiuid fiuxes 
are so much lower than in (a) and (c), that they are almost invisible, (e) Aquifer consolidation caused by the increase 
in effective stress following the tunnel construction (Weibull distribution). 



the aquifer consolidation increases because of the to- 
tal depressurization of the system (Figure jTjj). Over- 
all, the transient state is relatively fast because there 
is no release of water from the rock matrix, assumed 
impervious. 



Conclusions 

A model function relating fracture permeability to 
effective stress is derived from Hooke's law of elastic- 
ity and from the statistical distribution of asperity 
lengths. This model function is then implemented 
in the tensor form of Darcy's law, and its effects are 
assessed in simulations. Taking into account the sen- 
sitivity of permeabilities to depth and water pressure, 
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Figure 6: (a) Discharge rates, vertical settlement and (b) transit time as a function of coefficient n, and comparison 
with the classical approach for steady state flow before and after the tunnel construction. For transit time the particle 
is released at coordinates x — and z — 2300, and exits: (1) at the bottom of the valley (without tunnel); (2) at the 
tunnel (particle tracks are shown in Figure [s] for the classical approach and for the Weibull distribution) . 




Figure 7: Evolution of the discharge rate drained by the tunnel (for the classical approach and the linear increasing 
distribution, solid lines), and of the maximum vertical settlement (for the linear increasing distribution, dashed line) as 
a function of time for (a) constant hydraulic head at the domain surface; (b) no-flow condition at the domain surface. 



this non-linear approach gains in phenomenology and 
is closer to physical reality, compared to classical 
approaches that neglect pressure-dependent perme- 
ability and porosity fields. From a general point of 
view, numerical simulations of deep tunnels consid- 
ering the decrease in permeability with increasing ef- 
fective stress generate lower discharge rates. This 
observation goes in the line of safety in terms of 
problems due to the presence of water in the under- 
ground structure. In case of strong decrease in water 
pressures a non- negligible consolidation occurs, even 



when flow is assumed in the fractures only. This can 
produce foundation instabilities of structures located 
at the surface, especially in the case of differential 
consolidation due to aquifer heterogeneity, and in the 
presence of heavy structures such as dams (see [E]). 
The changes in fracture permeability, porosity and 
specific storage in response to changes in effective 
stresses depend on, among other factors, the statis- 
tical distribution of asperity lengths, indicating the 
relative ratio of large to small asperities. For exam- 
ple, in the case of increasing effective stresses, a frac- 
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ture characterized by a high ratio of large to small 
asperities (small n), will have a lower change in per- 
meability than a fracture with a small ratio (high 
n). Overall, the limitation of the classical method 
is that it cannot compute aquifer consolidation, be- 
cause no change in fracture porosity or permeability 
with pressure head variation is accounted for. On 
the contrary, with the proposed approach, a pressure 
head variation causing fracture porosity to change 
can be directly translated into aquifer consolidation 
(decreasing pressure) or expansion (increasing pres- 
sure) . 
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Appendix I: Notation 



a [m] 


fracture aperture (asperity length under compression) 


ao [m] 


maximum fracture aperture (original length of the longest asperities 


A [m^] 


area 


E [Pa] 


fractured rock elastic modulus 


d [m] 


distance 


D{z) [-] 


statistical distribution of the asperities length 


p r 1 1 

/ k ^] 


frequency of the fracture family 


F [AT] 


force 


g [m s~^] 


gravitational acceleration 


h [m] 


pressure head 


H [m] 


hydraulic head 


7 r On 

A; [m^J 


intrinsic permeability 


- r Oi 


maximum intrinsic permeability 




hydraulic conductivity 


T 7- r 1 1 


maximum hydraulic conductivity 


m [— ] 


number of fracture families 


n [-] 


coefficient of asperities length statistical distribution 


'^X; ) [ ] 


components of the unit normal vector 




number of compressed asperities 


^/ [-] 


number of fractures 


[-] 


total number of asperities 


p [Pa] 


water pressure 


s [m^] 


average asperity section 




specific storage coefficient 


[m-^\ 


maximum specific storage coefficient 


t [s] 


time 


T [m] 


ground settlement or expansion 


z [m] 


1: asperity's original length, 2: elevation head 


Z [m] 


depth 


a[-] 


Biot- Willis coefficient 


r] [m ^] 


asperity areal density 


A [-] 


ratio of horizontal to vertical stress 


[kg m s \ 


water viscosity 


Pr [kg m -"l 


rock mass density 


Pw [kg m 


water density 


a [Pa] 


normal stress 


a' [Pa] 


normal effective stress 


cTo [Pa] 


fracture closure stress 


a', [Pa] 


fracture closure effective stress 


0H 


porosity 


•/-o [-] 


maximum porosity 


A<p [-] 


porosity variation 
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Appendix II: Used values 





Fig. 3a 


Fig. 3b 


Illustration 
Rock 1 


Simulation 
Rock 2 


Rock 3 


m [-] 


1 


1 


3 


3 


2 


[mm] 




0.1 


0.5 


0.5 


0.5 


[mm] 






1.2 


1.92 


1.2 


aog [mm] 






1.0 


1.63 




fi [m-^] 




1 


5.44 


3.27 


0.1 








0.71 


0.71 


0.5 


h [m-^] 






1.00 


0.01 




Kmax [m S~^] 






1.56-10-3 


3.77- 10-3 


6.31 • lO-'' 


Kmin [m 






5.47 • lO-'' 


2.1 • 10-^ 


6.1 • 10-6 


9[°] 






27 


36 


34 




1.72 • 10"^'' 


9.24-10-^'^ 








a'o [MPa] 


350 


495 


350 


300 


325 


n[-] 


11 


2.5 


variable 


variable 


variable 


pr [kg m ^] 




2400 


2800 


2200 


2500 


f^xi 1T"yi "^z [ ] 




[1,0,0] 


[1,0,0] 


[1,0,0] 


[1,0,0] 








[-0.555,0,0.832] 


[-0.555,0,0.832] 


[-0.555,0,0.832] 








[0,0,1] 


[0,0,1] 




AH 




0.41 


1.5 


1.5 


1.5 


Z [m\ 




15 
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